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THE  MOTION  OF  A  SPAR  BUOY  IN  SWELL 
PART  II 

Philip  Rudnick 


Equations  of  motion  are  given  below  which  take  account  of  all 
the  pressure  forces  with  high  accuracy,  but  completely  neglect  drag 
forces.  The  results  of  the  previous  simplified  study  of  the  vertical 
motion  are  confirmed;  the  principal  new  result  is  the  amount  of  the 
horizontal  force  and  its  torque.  Nothing  large  or  important  is  seen 
in  the  interaction  terms. 


The  primary  dependent  variables  are  5  ,  n ,  and  8.  The  first 
two  are  displacements  of  the  center  of  mass  from  its  equilibrium 
position  in  still  water,  respectively  on  a  horizontal  line  along  the 
direction  of  wave  propagation, and  vertically.  Positive  directions 
are  respectively  with  the  wave  travel,  and  upward.  0  is  an  angle  of 
rotation  about  a  horizontal  axis  parallel  to  the  wave  fronts.  Positive 


sense  for  0  is  from  the  positive  5  direction  to  the  positive  n  direction. 


Notation:  Upper  and  lower 
cross-sections  are  rA  and  A, 

Axial  positions  are  measured  from 
point  P  which  is  at  the  waterline 
when  the  buoy  is  at  rest  in  still 
water.  Meanings  offc^  and£2 
are  indicated  in  the  diagram; 

Ji,  and  %  give  respective  positions 
01  cental-  of  buoyancy  (for  immer¬ 
sion  in  still  water  to  point  P)  and 
center  of  mass, 

1 

Mass  of  the  buoy  is  expressed 
as  P  AL,  where  p  is  density  of  the 
displaced  sea  water.  We  have 
L  =  *2  -  (1  -  r)A  .  R  is  the  radius 
of  gyration  about  the  center  of  mass 


\ 


A  Area 
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Pf  is  the  instantaneous  intersection  of  the  buoy  axis  with  the 
undisturbed  sea  surface  and  P"  its  intersection  with  the  actual  sea 
surface.  P'  and  P"  move  with  respect  to  the  buoy,  while  P  is  fixed. 

6  andc  measure,  the  displacements,  positive  when  upward,  PP*  and 
P’P"  respectively.  0  is  the  inclination  of  the  sea  surface  at  P”  from 
the  horizontal,  with  the  same  positive  sense  as  0, 

Description  of  wave  motion:  Eulerian  positions  x  and  y  will  be 
referred  to  the  same  coordinate  directions  as  5  and  n  ,  respectively. 

The  undisturbed  sea  surface  is  y  =  0.  The  waves  will  be  assumed 
to  be  unidirectional,  in  the  positive  direction,  with  surface  amplitude 
a  and  single  radian  frequency  oi  .  The  corresponding  radian  wave 
number  k  is  given  by  u2  =  gk  with  g  the  acceleration  of  gravity. 

Particle  displacements  x*  andy*  are  then  given  as  functions  of 
position  (x,  y)  and  time  t  by 

ky 

x*  =  ae  3  cos  (kx  -  cot) 

y'  =  ae^  sin  (kx  -  cot)  (1) 

The  usual  expression  taken  for  the  pressure  p  is 

ky 

p  =  pg  [  -y  +  ae  3  sin  (kx  -  tot)]  (2) 

(1)  and  (2)  are  linear  approximations;  the  limitation  in  accuracy  can 
be  seen  by  noting  that  from  (1)  the  height  of  the  sea  surface  is 

y^  =  a  sin  (kx  -  u  t) 

whereas,  from  (2)  the  locus  of  vanishing  pressure  is  the  solution  of 

ky 

y  =  ae  3  sin  (kx  -  w  t) 

lest 

These  agree  only  if  a  factor  as  large  as  e  is  ignored  in  the  exponential. 
Hence  (2)  may  be  and  was  replaced  by 

k(y  —  y  * ) 

p  =  P  gC  -y  +  ae  3  3 o'  sin  (kx  -  w  t)  ]  (3) 

which  vanishes  accurately  at  the  sea  surface.  Components  of  the 
pressure  gradient  were  taken  as 
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±2. 
3  X 


k(y  -  y* ) 

gaku  7  7o  cos  (kx  -  w  t) 


3 -ft  =  p  gC  -1  +  ake^y  V  sin  (kx  -  o>  t)  ] 

ay 

The  value  of  the  phase  kx  -  to  t  at  point  P"  is  abbreviated  to  u  .  We 
also  have 


(4) 


6  =  Z  g  (sec  0  -  1)  -  n  sec  0 
e  =  a  sec  0  cos  a 
0  -  ka  cos  a 

General  description  of  calculation:  The  lateral  pressure  force  on 
a  volume  element  of  the  buoy  whose  base  area  is  A  and  whose  axial 
extent  is- d£  was  taken  to  be  the  product  of  the  volume  Adi  and  the 
component  of  the  pressure  gradient  (4)  on  the  buoy  axis  and  normal  to  it. 
Total  lateral  force  or  its  moment  was  calculated  by  integration  of 
such  elements  between  P”  and  the  bottom  of  the  buoy.  The  fact  that 
the  sea  surface  was  not  necessarily  normal  to  the  buoy  axis  at  P"  was 
considered  to  require  no  correction  to  die  lateral  force/  and  only 
a  small  one  to  the  torque.  The  reasoning  is  indicated  in  connection 
with  the  diagram  shown  below.  The  integration  over  all  sections 
up  to  P"  includes  the  set  of  (negative  pressure)  forces  which  does 
not  actually  act,  and  fails  to  include  the  set  F2  which  does.  Fj  and  F2 
are  taken  to  have  equal  resultants,  so  that  the  only  effect  of  replacing 
one  by  the  other  is  a  small  torque  depending  on  their  difference  in 
axial  positions.  / 
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Results:  The  differential  equations  of  motion  can  be  written 


~2  5  =  “g<a/L)h 


(d2n  /dt2)+u)2n=w  2  nJl  -  (ga/L)  v 


(d2  0/dt2)  +  to  2  sin  9  =  (-gr/LR2)  (sin  0)  ^  +  (ga/LR2)  T2 
+  (gr/LR2)  T3 


in  which 


h  =  r  cos  (a+  0)  +  (1  -  r)  e  ^  COS  ®  +  <5+  e  )cos  (  a+  0+  ksin  0 

(4  +5+e  }  )  -e”kCOS  0(*2+6  +E  ^cos  (a+  0  +  k  sin  8{f*+6+e  }) 
*  [  r  +  (1  -  r)  e  1  -e  ?  ]  .cob  a  (8) 


a)2  =  (gr/L)  w2  =  g  (£g  -  y  /  R2 


nl  =  "  n)  <sec  0  -  1)  -a  tan  0  cos  £  +  0) 

O 

-  (1/2)  4  02  -  a  0  cos  a 
g 


=  (1  -  t)  e  ^  C°S  ®  ^1  +fi  e  ^sin  (a+  0+ksin0{£1  +  6+e  }) 
-e  k  cos  0  (  +  <$  +e  ^sin(a+0  +  ksin0{£2-hS+e}) 

[  (1  -  r)  e  ^*1  -e  ^  2  ]  gin  a  (11) 


c-p  ,  .  1  v  12 

T1  =  fcg  +  2  5)5_2  ^ 

*  -4  n+^C(40)2  +n2  -(asino)2]  (12) 

g  ^ 
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T9  =  r[(£r+5+e)  sec  0  -  1/K  jcos  (<*+  0)+(l -r)e"k  COS0(*  1  +  *  +  6  } 

[  (  t  r  a  cos  (  ot  +  k  sin  9U^  +  5  +  e}  )  -  1/k  cos  (  a+  0  +  k  sin  0  +  e}  ) . 

+  e  k  C0S  ®  ^2  +  <S+t  ^  [  (£  -  £  )  cos  (a+  k  sin  G{£  0  +  5  +  e}  ) 

2  g  2 

+  1/k  cos  (  a  +  0  +  k  sin  8  { £  ^  +6  +e  }  )] 

-ko  -k  £ 

=•  [  r(n  -  l/k>+  (1  -  r)  c  "  1  ( ,  -  t  -  l/k)+e  2  ( i  - 1  +  l/kjjcos  „ 

O  O  o 

-(£„-  £n)h  (13) 

g  0 

T  --  (1/8  "  )  r  A  sin  (0  -  0)  tan^  (0-0)  =  0  (14) 

o 


In  all  of  the  expressions  (8)  through  (14)  there  is  given  first  an 
instantaneous  value  as  accurately  as  it  was  calculated.  This  is  followed 
by  a  suitable  simplification  in  which  some  small  terms  are  dropped, 
and  others  altered  to  be  applicable  to  a  whole  cycle  of  motion,  assumed 
periodic  with  fundamental  frequency  w  .  In  the  simplified  expressions  a 
must  be  understood  to  refer  to  the  mean  position  of  PM  throughout  a  cycle. 

It  is  evident  that  to  a  good  approximation  the  three  coordinates  may 
be  considered  separately.  By  (5),£  may  be  considered  as  a  forced 
oscillation  of  a  free  particle. 

n  may  be  considered  as  the  motion  of  an  oscillator  of  natural 
frequency  ^  under  a  periodic  impressed  force  represented  by  the 
term  containing  v  in  (6).  The  terms  in  n  ^  are  quadratic  and  hence  tend 
to  contribute  small  constant  displacements  and  second  harmonic 
distortion. 

Similarly  the  motion  in  0  is  largely  governed  by  the  natural 
frequency  wQ  and  the  periodic  impressed  torque  represented  by  the  term 
containing  in  (7).  The  term  involving  T^  arises  from  lack  of  per¬ 
pendicularity  between  die  buoy  axis  and  the  sea  surface,  and  is  quite 
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small.  A  positive  constant  part  of  Tj  represents  increased  restoring 
torque  and  a  consequent  rise  in  natural  frequency.  The  periodic  parts 
of  Tj  contribute  constant  displacement,  or  second  or  third  harmonic 
distortion.  These  perturbations  will  all  generally  be  small. 

The  first-order  steady  state  motions  are 


=  (a/kL)  h 

2  2 

(15) 

n 

=  [  ga/L  (  w  -  o>v)  ]  v 

=  av  /  (kL  -  r) 

(16) 

e 

=  { (£  -  £  Q)  /  [  R2  -  ~  (  £  - 
q  W  k  g 

tj]  }  (a/kL)h 
b  (17) 

The  horizontal  amplitude  approaches  a  for  small  k  fand  to ),  and 
is  considerably  less  than  a  at  higher  frequencies.  The  vertical  ampli¬ 
tude  has  a  resonant  peak  when  kL  =  r  and  a  zero  when  1  -  r  =£“k(  2  ^ 
The  zero  will  always  occur  at  a  higher  frequency  than  the  resonance. 

The  vertical  amplitude  also  approaches  a  at  low  frequency,  and  vanishes 
exponential) v  at  high  freau°ncy. 

The  distance  £g  defired  in  (13)  locates  the  effective  center  of 
horizontal  wave  force  relative  to  point  P.  Its  value  varies  from  zero 
at  very  high  frequencies  to  £  ^  at  very  low  frequencies.  The  relation 
(17)  shows  directly  the  relation  between  angular  displacement  and 
linear  displacement  at  die  center  of  mass.  The  center  of  oscillation  can 
be  varied  over  quite  wide  limits  by  choice  of  the  resonant  frequency  ui  q. 

Effects  of  drag:  The  velocities  indicated  above  for  the  buoy  are  usually 
in  phase  with  (but  sometimes  differ  by  13C°  chase  from)  the  water  wave 
velocities.  In  either  case  the  drag  forces  will  add  energy  to  the  motion 
where  the  water  velocity  is  in  phase  with,  and  greater  than,  the  buoy 
velocity;  otherwise  damping  will  cccur. 

The  following  order-of-magnitude  estimates  of  the  ratio  of 
drag  to  pressure -plus -gravity  force  have  been  made  on  the  assumption 
that  the  buoy  is  stationary.  For  the  horizontal  motion,  assuming 
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unity  drag  coefficient  for  a  cylindrical  form  against  transverse  flow, 
the  ratio  at  any  frequency  is  of  the  order  a/D.  where  D  is  the  diameter 
of  the  cylinder.  For  end-drag  due  to  vertical  flow,  at  any  frequency 
the  ratio  is  of  the  order  of  ka,  again  assuming  unity  drag  coefficient. 
For  skin-  drag  arising  from  vertical  flow,  the  ratio  at  high  frequency 
is  of  the  order  (ac)  /  (Dv)  where  c  is  the  skin-drag  coefficient,  with 
a  possible  value  in  the  range  0. 002  to  0. 003.  At  low  frequencies,  the 
ratio  is  less  than  that  given  by  this  expression.  This  result  indicates 
a  lower  limit  for  die  vertical  pressure  excitation,  below  which  it  will 
be  superseded  by  skin  friction  excitation. 
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